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O ' Abstract 

We treat a system (a molecule or a solid) in which electrons are coupled 
linearly to any number and type of harmonic oscillators and which is fur- 
-^ , ther subject to external forces of arbitrary symmetry. With the treatment 

[^^ I restricted to the lowest pair of electronic states, approximate "vibronic" 

^O ■ (vibration-electronic) ground state wave functions are constructed having 

the form of simple, closed expressions. The basis of the method is to re- 
(^ I gard electronic density operators as classical variables. It extends an earlier 

"guessed solution", devised for the dynamical Jahn- Teller effect in cubic 
symmetry, to situations having lower (e.g., dihedral) symmetry or with- 
out any symmetry at all. While the proposed solution is expected to be 
'^ ' quite close to the exact one, its formal simplicity allows straightforward 

^ ■ calculations of several interesting quantities, like energies and vibronic re- 

Q , duction (or Ham) factors. We calculate for dihedral symmetry two different 

g-factors ("Qz" and "g^") and a p-factor. In simplified situations we obtain 

p = qz + Qx-l- 
The formalism enables quantitative estimates to be made for the dy- 
d I namical narrowing of hyperfine lines in the observed ESR spectrum of the 

dihedral cyclobutane radical cation. 



PACS: 71.70.Ej, 31.30.Gs, 76.30. 



1 Historical Background and Aims 

For the so called E ® e Jahn- Teller case (involving an electron-nuclei sys- 
tem, in which a doubly degenerate electronic state is coupled to a doubly 
degenerate nuclear displacement mode) the wave function was fully obtained 
as long ago as 1957 [1, 2]. The physical object of reference is commonly a 
molecule of some high symmetry (say, one belonging to the cubic group, like 
O), or a localized impurity in a solid. The solution (or set of solutions) to 
this "dynamic Jahn- Teller effect" (DJTE) are the vibronic states. Though 
this has received, as just noted, a full treatment early on, subsequent efforts 
to give simple approximate treatments or to provide additional insight into 
the dynamic problem have been numerous. Descriptions of some of the early 
works are found in two books [3, 4]. Notable are the treatments in [5]-[9]; the 
most recent publication known to us and involving a variational approach 
to this problem is in [10]. 

To lead us into the present work we recall a "guessed solution" for the 
ground state of the linear Jahn- Teller effect, suggested by one of the present 
authors and collaborators, which is transparent, intuitively simple and al- 
gebraically easily manageable. This proposal was originally worked out for 
a molecule of cubic symmetry which had a single set of interacting normal 
modes [11, 12, 3]. Though not variationally obtained, the "guessed solu- 
tion" was found to have energies that are considerably closer to the exact, 
computed energies of [2] than other approximate solutions with which it was 
compared. This comparison is seen in Fig. 2 of [13]. Later treatments did 
not test their methods by comparison with the "guessed solution", though 
a critical review can be found in section 4.5.3 of [4]. 

The present work is an extension of the earlier approach to a substan- 
tially broader and harder problem, namely to a pair of electronic states in 
unrestricted symmetry and subject to interaction with an arbitrary number 
of nuclear displacement modes, but only in a linear manner. The subject of 
two-state interacting with bosons (which may either be phonons or photons) 
has had a very extensive literature. The spin-boson Hamiltonian that forms 
the starting point of [14] is a special case of the Hamiltonian introduced in 
this paper. Likewise, several books contain accounts of the related Jaynes- 
Cummings method [15, 16]. The present work also belongs to this field, but 
is restricted to a pair of ground level states. Even with this restriction, the 
closed solution that we present here can find its uses in treating the energy 
dissipation of a spin system [14]. 

The handling of external perturbation after taking care of the electron- 



nuclear interaction is a potential tool to tackle Berry phases in open systems 
[17, 18]. We would also recall a recent work on the Jahn- Teller effect in 
lower than cubic symmetry, which is less general than the present one, but 
has permitted us to check some of our results numerically [19]. The reduced 
symmetry case (named "the elliptic form" to differentiate it from the circular 
energy trough in £? e) was studied previously in [20] . We calculate (for the 
first time, to our knowledge) the experimentally important reduction factors 
for the low symmetry case (section 4.1), having pointed out (at the end of 
section 3.4) that, when the electron-nuclear coupling is strong, one meets 
broken symmetry instabilities. 

The formalism, initially formulated in very general terms, is gradually 
shifted to more specific situations, such as systems of cubic and of lower 
(e.g., dihedral) symmetries, and to systems with two (rather than an arbi- 
trary number of) vibrational modes and, ultimately, to a specific molecular 
system. In this last, the formalism and the numerical results for the reduc- 
tion factors lead to quantitative conclusions for the dynamical narrowing of 
hyperfine lines in the observed ESR spectrum of the dihedral cyclobutane 
radical cation. This is the subject of section 5.2. 

2 A General Hamiltonian 

We now write down a Hamiltonian for a pair of (diabatic, or nuclear 
coordinate-independent) electronic states, denoted by the symbols ^6* and 
Ce). The states are understood to be functions of any number of electronic 
coordinates, e.g, all the electrons in an atom or a molecule, but this func- 
tional dependence is absent in the formalism as long as the behavior of the 
doublet states alone is under consideration. A two-state situation can come 
about for an atom that is placed in a strongly coupled environment, such 
that this separates the doublet from the rest of the electronic manifold, or 
for a molecule in which the internal, intramolecular forces achieve the same 
effect. The two states need not be degenerate but, in order that it should be 
legitimate to consider them separately, they must be, in some sense, isolated 
from the rest of the electronic states (e.g, either by symmetry consideration 
or by a large energy gap). If so, then the "two level-the rest" matrix ele- 
ments of all interactions can be neglected to some approximation. This is 
the physical setting for the formalism that follows. It leads naturally to the 



representation of the electronic states as the column vectors 



The electronic states are coupled to any number of nuclear displacements 
coordinates g„ (n = 1, ...N). We assume that these are organized into a set 
of normal modes brought to a standard form (i.e., having the same effective 
mass) and restrict the coupling to be of no higher order than linear in the 
displacement coordinates. Thus, one has for the displacement coordinates 
the following harmonic oscillator Hamiltonian: 
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where hujn are the quanta of vibrational energies. In the two-state repre- 
sentation the nuclear Hamiltonian is written as a scalar or, equivalently, as 
Hnuc times the 2x2 unit matrix /. 

The remainder of the Hermitian representation matrices for the two level 
system are the familiar Pauli-matrices 



a = ia,,ay,a,) = {[ , r, ] , [ , n ' n _i ) (3) 




In terms of these we can write out a general linear form of interaction be- 
tween the electronic motion and the (real) nuclear coordinates (in the ab- 
sence of any molecular symmetry) , as 

Hel-nuc = - X! ~^i'^nqn(yz + Klnf^x) (4) 

n=l ^ 

in which (the dimensionless) a^ and 6„ express the strength of interaction 
between the electrons and the nuclear motion in the n- mode. (A detailed 
discussion of the linear many mode interaction in a symmetrical setting is 
found in section 3.5.3 of [4]. In equation (4) frequency changes between the 
two states are ignored to be consistent with a purely linear coupling.) 

(Ty is absent in the above interaction Hamiltonian, as also in several pre- 
vious works [l]-[8]. When the states of the two level system are orbital states 
then, for all molecular point groups considered in this work, the symmetric 
product of the state-representations does not contain the representation of 



the cjy-niatrix. When the two states are a Kramers-doublet, the situation 
becomes more complex, since in a low order perturbation the coefficients a„ 
and bn vanish, unless some further effects (like crystal field, spin-orbit inter- 
action, external magnetic fields H, spin-spin coupling) are included in the 
perturbational calculation of these coefficients. Working out the spin-lattice 
coupling for a Kramers' doublet (on a six coordinated Cu^~^) Stoneham gave 
symmetry arguments (in the last equation of [23]) to show that, for both 
a„ and 6„ to be non-zero, both H^ and Hz need to be non-vanishing (while 
Hy = 0). This corresponds to the form shown in the above equation (where 
the coefficients would be magnetic field dependent). If, on the other hand, 
Hy is also non- vanishing, then for a Kramers' doublet a term with ay will 
also be present. Since this term is absent in orbitally two-state system, we 
do not complicate the formalism by adding the ay term. 

One also has to consider the electron being acted upon by external fields. 
(The interaction of the external field on the nucleus is supposed to be con- 
tained in the potential of the nuclear coordinates.) In the preceding, vector 
representation of the two states, any interaction Hamiltonian (that expresses 
the coupling between the electron and any external field) must have the form 

Hel-f = 1^ • <? (5) 

with the representative of the fields il = {i^x,^y,^z) inside the two level 
system being constant (independent of the values of the electronic or of the 
nuclear variables), this being the most general form of expression for the 
system. In section 4, which discusses the effect of external forces, we give 
examples for the interaction. 

Any difference between the two-state energies can be considered to be 
part of Qz so that, until we come to the subject of the external fields in 
section 4, the states can be considered as a pair of degenerate doublets. 
However, the rest of Qz (as well as Qx and fi^) comes from externally applied 
sources. 

The total Hamiltonian Htot is the sum of the previous Hamiltonians 

Htot = {Eo + Hnuc)I + Hel-nuc + Hel-f (6) 

to which has been added a scalar term with Eq representing the mean energy 
of the non-interacting states. (The spin-boson Hamiltonian which forms the 
basis of [14] is obtained from equation (4) and equation (5) upon putting 
a„ / 0, fe„ = 0, [7^ / / Qz,% = 0.) 



The treatment of H^i-f will be postponed to later. In its absence, we 
have a pure "vibronic (=vibrational-electronic)" situation, which we now 
treat. 

3 Vibronic doublet 

The Hamiltonian 

^v — J^nuc ~r -tlel—nuc ( ' ) 

involving (partially) the electronic and nuclear degrees of freedom will be the 
subject of our investigation in this section. We first show that solutions of 
the partial Hamiltonian form a degenerate doublet (with the understanding 
that the energy difference between the two states is shifted to the external 
field part). 

3.1 State degeneracy 

The following non-identity transformation leaves the above Hamiltonian in- 
variant: 

'r = [Ae][r,] (8) 

where Ae stands for the following simultaneous changes in the electronic 
for all mode coordinates, namely 



states I ri 1 ^ I „■ 1 ^^'^ I 1 1 ^ I n 1 ^^'^ ^q i^ ^^^ parity operator 



Pqqn = -qn u = I, ...N (9) 

The transformation Ae can be achieved by the unitary matrix ay. But, 
clearly, 

r^ = 1 (10) 

so that, simultaneously with any eigenvalue of the Hamiltonian H^, T 
has an eigenvalue -|-1 corresponding to a state unchanged under the T- 
transformation and another eigenvalue —1, that corresponds to (another, 
different) state that changes sign under this transformation. The conclusion 
is that the vibronic Hamiltonian has doubly degenerate eigenstates. This 
degeneracy is lifted, when field interaction term H^i^j is inserted, as we 
shall see. 



3.2 The quasi-classical ground state 

Based mainly on the numerical agreement of the energies under O symmetry, 
noted in the opening section, we extend here the method of [3], [12] to the 
general case under study, that is, we propose the following form for the 
ground state wave-function 

1 ^ h 

^ = AA'exp{-i ;^(g„/ - ^a. - |a.)n (11) 

n=l 

where N' is a normalizing factor. The wave function generator \E' is a ma- 
trix (or operator). It possesses the full {Ai) symmetry of the Hamiltonian. 
Therefore operating with ^ on an electronic component with some symme- 
try, will generate a state with the symmetry of the component. In order 
to get the probability amplitudes in the two electronic states of equation 
(1) , one has to left-operate with ^ on the basic vectors or on any linear 
combination of them, as will be shortly described. The prescription (and the 
underlying rationale) for the proposed construction is to regard the Pauli 
matrices (which are the electronic density operators) as c-numbers. Having 
done this, we write down the ground state wave function in the form of a 
set of displaced independent vibrational coordinates. The initial handling of 
(quantum mechanical) matrices in the manner of c-numbers has suggested 
naming the method "quasi-classical". However, the modes are now no longer 
independent: thus the moments (e.g., the expectation value or the spread) 
of any mode depends on the coupling constants of the other modes. 

The mathematical meaning of the exponential form in equation (11) is 
that one has to expand the exponential in a power-series of the exponent. 
Remarkable in the posited ^ is that the individual frequencies do not appear 
in it (just as they do not in the wave function of a set of uncoupled oscillators, 
when expressed in a standard form). Of course, the energy expectation value 
depends on the frequencies, since the Hamiltonian does and so do, implicitly, 
the dimensionless coupling constants a„ and 6„. 

The success of the method hinges on the fact that it is possible to sum the 
power series exactly, in spite of the non-commuting terms in the exponent. 
This is made possible by the property of the two dimensional spin (Pauli) 
matrices, that 

aiUj + ajGi = 25ij (12) 

where 5ij is the Kronecker delta. Using this property, one readily obtains 



the simplified expression 

^ = N[I cosh £M + Ell(^.angn + a.bngn) e(^ 

with A^ being another normaHzing factor. In the argument of the hyperbolic 
functions one has 



e(f) 



\ 



N N 

[{Y.^nan? + {Y.^nbn?\ (14) 



It must stressed again that the resulting quasi-classical wave-function 
is only an approximation, whose accuracy depends on how well the Pauli- 
matrices can be approximated by c-numbers. This will be the case when 
e.g., one of the potential wells is deep, since then the energy of the state will 
be dominated by the electron occupancy near the minimum. On the other 
hand, when the frequencies of the different oscillators differ markedly, the 
quasi-classical approximation could be in error, since then, e.g., the positions 
of the saddle points in the potential may not coincide with the maxima in 
the overlap of wave-functions coming from different potential wells. A scale 
transformation applied to each well, in the form proposed in [10, 21, 22] for 
cases of higher degeneracies than two and modes of various dimensionalities, 
could lead to improvements in the wave-function, but requires a formalism 
that is more complex than the one advocated here. 

3.3 Some elementary symmetry considerations 

As already emphasized, there need not be any relation (symmetry-based or 
otherwise) between the two states and among the nuclear coordinates for 
the foregoing formalism to hold. However, if the system has some symme- 
try properties (or we choose to relate it to a symmetric framework) things 
become at the same time clearer, more systematic and more familiar. 

We therefore formulate the foregoing in a symmetry-setting and employ 
implicitly the theory of point molecular or point-groups. 

3.3.1 Cubic symmetry 

In a system which nominally belongs to a cubic symmetry group (like O), the 
two electronic states could belong to a doubly degenerate E representation, 
whose two components are designated in [24] as (0, e). This designation was 



already used by us above in equation (1) . We have chosen these symbols 
in preference to others, such as those in [25], because we shall later use the 
so-called W-coefficients and these were tabulated for the symmetry groups 
of interest in [24] using the present symbolism. In the expression, equation 
(4) , for the coupling to the nuclear motion coordinates g„, the coefficients 
a„ are non-zero for modes belonging to the ^-representatives of a two-fold 
e-mode and 6„ are non-zero for modes belonging to the e-component of an 
e-mode, with the two coefficients being numerically pairwise equal. This is, 
of course, the multi-mode £' (gi (ei + 62 + ... + e^/2) Jahn- Teller situation, 
described in detail in [3], [4]. 

^ in equation (11) generates vibronic wave function in the following way. 

When it is let to operate on , one obtains the 9- component of the 

ground state vibronic doublet (in the present quasi-classical approximation) ; 

if it is let to operate on , one obtains the e-component of the same. 

For future use we write these vibronic wave-function components in a curly 
ket form, as 

^(l]=\e}, ^(l]=\e} (15) 

If, instead, one operates with ^ on the two orthogonal linear combina- 
tions J- 1 1 ' ^'^^ reaches a pair of other vibronic states, preferentially 

localized in a different part of the coordinate space than the ones in equa- 
tion (15) . (These vibronic states have properties similar to polaronic states, 
which term is in use for a single electronic state.) 

If, alternatively, one operates on the following complex combinations of 

the electronic states, r- • , one obtains vibronic states that have 

values of +1 and —1 of a "composite" vibrational angular momenta, this 
being defined in terms of the composite vibrational angle variable, given by 

q? = arctan j^ (16) 

Clearly, although the vibrational modes were originally independent, the 



resulting angular variable is not the sum of the mode angular variables, like 



JV 

- N 



<& = ^ ^ arctan ^^ (not true) (17) 



m=i™/=4+i ^"^ 



2 



This is, of course, due to the coupling of the modes to the electronic degree 
of freedom. 

A pair of states in any of the combinations are energy-degenerate and 
mutually orthogonal. 

3.3.2 Dihedral symmetry 

In a lower symmetry situation, like D2d, the doublet state could belong 
to the doubly degenerate i?-representation, with components designated by 
{0,e), as before. The a„ coefficients are non-zero for modes possessing 62 
symmetry and the 6„ coefficients belong to 61 symmetry types. However, 
this time (unlike for 0-symmetry) there are no symmetry-based relationships 
between a's and 6's. The operator ^ still possesses full (^1) symmetry. This 
can best be seen in the form given in equation (13) . Here the first term 
is clearly invariant under all group operations (since the squares of these 
are the identity operator) and the second term is also invariant because it 
has the symmetry of the Hamiltonian. Therefore operating with ^ on an 
electronic component with some symmetry, will again give a state with the 
symmetry of the component. 

3.4 Energies and other expectation values in dihedral sym- 
metry. 

The advantage of the quasi-classical form is that the expectation values of 
c-numbers or of operators can be calculated using the explicit form of the 
states just given. (A recent calculation of the expectation value of an angular 
momentum operator is in [26].) As already noted in the opening section, 
these expectation values have proven to be very accurate for the single mode 
case in O symmetry, where comparison was made with the exact, computed 
results that were available. One can carry out similar calculations for several 
modes in any dihedral group, like one of D4 symmetry (as well as in other 
dihedral groups, like D^h, D2d, etc). Actual results will not be given here 
for several modes, since these will depend on details of the system. (Some 
relevant molecular systems will be considered in section 5.) 
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A sporadic comparison has been made with results of a recent paper 
[19]. This paper showed graphs of eigenenergies of the (degenerate) vibronic 
ground computed exactly (by a numerical method), as well as with several 
approximation schemes. We compare the energy expectation values com- 
puted within our quasi-classical approach with theirs, at values of parame- 
ters for which the discrepancies between the exact and approximate values 
appear to be largest. This is the region where neither perturbation theory 
(weak coupling), nor asymptotic formula (very strong coupling) holds. Re- 
lating to figures 3(a) and (b) in [19] and to parameter values a = 2, 6 = 1.5 
(equal to fi = 0.5, x = 0.75 in the symbols of [19]), their exact eigenenergy 
IS E = —0.16725. We compare this to —0.15075 obtained by our method, 
which is higher (as it should be for a non-exact expectation value) by 0.0165. 
This discrepancy is worse than the best approximation in [19], rather better 
than the next best (obtained variationally) and considerably better than 
three others. Testing additionally our quasi-classical approximation against 
the exact results exhibited in the figures 4 (a) and (b) of the above ref- 
erence, for a = 2\/2,b = 1.5\/2 {/j, = 1.0, x = 0.75), the exact value is 
E = —1.17295, with which we can compare our value of —1.13370, or a 
discrepancy of 0.03925. This is somewhat worse than the first and second 
best approximations of [19]. 

In summary, it seems that the intuitive, quasi-classical method can an- 
swer most practical needs for energy values. One would expect at least semi- 
quantitative guidance for other numerical quantities, when derived from the 
quasi-classical model. The next section contains some such quantities. 

It is also of interest to consider the limiting case of very strong electron- 
vibrational coupling. This comes about when at least some of the coupling 
strengths a„ and 6„ in equation (4) are numerically much larger than unity. 
(The opposite extreme of zero coupling has for its "vibronic" ground states 
the product of gaussian states multiplying some linear combination of the 
electronic basic states.) 

In the strong coupling limit the vibronic ground state is stabilized by the 
coupling and again takes the form of a product of vibrational and electronic 
factors. The form of the latter depends critically on the ratio of two stabi- 
lization energies. (This result goes back to Opik and Pryce's classic paper 
[29]. Modifications due to higher order coupling are treated in [30].) The 
stabilization energies (analogous to the static Jahn- Teller energy [3, 4]) are 

1 ^ 
AEb, = - ^ hiJnianf (18) 

n=l 
11 



1 ^ 

AEb, = ^J2 ^n{bn? (19) 

n=l 

where the subscripts 62 and 61 signify the types of the distortion mode in 
D2d symmetry. The ratio of the two 

enters now so that for R > \ the stabihzed electronic states are the {6, e) 
components, shown in equation (1) , and the stabihzation energy is A£^(,2) 

while for R < 1, the stabilized electronic states are i- , and the 

stabilization energy is AE'f,^. (i? = 1 is a coincidental case in dihedral 
symmetry. It is the normal case under cubic symmetry and results in a 
continuum of stable configurations, rather than a single point in the multi- 
dimensional configuration space for each state.) 

The transition between R < 1 and R > 1 represents a change from 
one broken symmetry type to another. The subject of broken symmetries 
is, of course, an important issue for macroscopic systems and also for the 
observed inhomogeneous state of the universe. In these cases each stabiliza- 
tion energy is large on a characteristic quantum scale (in a macroscopic or 
in an astronomical manner). Yet, it is a tiny difference between the sta- 
bilization energies that tips the balance between different types of broken 
symmetry. This means that if, in the neighborhood of a situation where the 
stabilization energies are the same, the coupling constants are made varied, 
then changes in the macroscopic symmetry and energy can come about by 
microscopic causes. (This is exemplified in our treatment of the cyclobu- 
tane radical cation in section 5.2.2.) The foregoing treatment is of course 
very approximate and does not take into account higher order couplings or 
thermal fluctuations; still, as a model it is instructive and can prove to be 
helpful for extensions to more realistic cases. 

4 External Fields 

The form of the interaction, shown in equation (5) , between the electronic 
part and an external field (represented by Q) is the consequence of the 
Hermitian nature of the Hamiltonian, the three matrices (cJ^, o^j/jCr^) being 
the only 2x2 matrices having this property (apart from the unit matrix, 
which only shifts both states by an equal, constant amount). 

12 



4.1 Examples for Vl 

The simplest example is the Zeeman effect on an electronic spin, for which 
O = PH (with P being the Bohr magneton). The effect of a magnetic field 
Hz acting within the two states {t2g,(,t2g,r^), which are split off from a 3d- 
state manifold by crystal fields of cubic and tetragonal symmetries, would 
be represented by 0,y =constant H^, where the constant includes /3 and a 
radial integral ([23], Table 1). A uniform stress of the type Tzz acting upon a 
doubly degenerate set under C2v, with z the fourfold axis, will be expressed 
by i^z- An applied variable electric field given by a potential V{r) = a;y/(|r|), 
when acting within an electronic pair having the real forms |1 >= 2;(7(|r|), 
|2 >= yh{\r\) will yield 

f^^ =< i|y(f)|2 > (21) 

with all other fi-components being zero. 

In a general way, for a perturbational Hamiltonian AH being an arbi- 
trary function of the coordinates, the fi- magnitudes and real and imaginary 
parts of the matrix elements inside the 1,2 manifold are connected by 

a, = i?e(AF)i,2, fty = -Im{AH)i,2, ^z = ^((AF)i,i - (AF)2,2) (22) 

Since the interaction term H(,i^f leads naturally to consideration of the 
vibronic, or Ham-, reduction factors [27, 28], we introduce these now. 

4.2 Reduction factors 

One starts with some coupling affecting the electrons. Within the two state 
manifold, this can be expressed in terms of matrix elements between these 
(adiabatic) states. The question is how do these matrix elements change, 
when the two states are no longer purely electronic, but rather vibronic 
(coupled electronic-vibrational) states? The answer, given in a context sim- 
ilar to the present one and originally due to Ham [27, 28], is a reduction 
in the strength of the original coupling by factors originally denoted by q 
and p (that are 1 in the absence of vibrational coupling and less than 1 
in their presence). These reduction factors have further been described in 
[3] and, at considerable length, in [4] where references to several literature 
sources can be found. Section 4.7 in the book [4] contains an analysis of the 
reduction factors in terms of W-coefficients, tables for which (in both cubic 
and dihedral groups) can be found in [24]. 
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In the lower symmetry situation under concern here (e.g., in dihedral 
symmetry) , when the coupling strengths in the diagonal position differ from 
those in the off-diagonal positions, there are three reduction factors. These 
are here named qx, p and Qz- The formal definitions are in terms of the 
scalar product of the vibronic states introduced in equation (15) , in the 
form 

qx = {9\cTx\e} = {e\ax\9} (23) 

qz = {eWM = -{^WM (24) 

p = i{e\ay\e} = -i{e\ay\9} (25) 

(On notation: The subscripts of the g's agree with those of the cr-matrices; 
the letter p has been retained in preference to a possible iqy for historic 
reasons. In terms of the more recent symbols K{a) used in, e.g., section 
4.7.1 of [4], where a is a representation of the group, one can identify: 

Qx = K{B2) , qz = K{Bi) , p = K{A2) , 1 = K{Ai) (26) 

The symbol K{a) will be used later in obtaining formal expressions for the 
reduction factors. In a higher symmetry situation (e.g. O) qx = qz = q, as 
shown on p. 43 of [3], where the ket in Eq. 3.38 should be corrected to be an e- 
type.) The results of our computations are shown for the simplified situation 
that there is a single 62 and a single 61 mode coupled to the electrons. When 
the coupling strength are of the same strength (a = ±5), our results are 
identical to those quoted in the literature for O-symmetry. In particular, 

qx = qz = q, p = 2q — 1 {in O symmetry) (27) 

but when the coupling strength are unequal a remarkable change occurs, as 
shown in the following figures (Fig. 1,2). 

These exhibit the three reduction factors as function of increasing 
strength of the dominant coupling and for two different values of the ratio 
b/a. It is apparent that while one (7-factor decreases monotonically to zero 
(!), the other drops only slightly below unity and does so for only a lim- 
ited range. (In cubic symmetry, also for linear coupling, the single g-factor 
decreases from unity to one-half.) The p-factor shows a regular behavior. 

The following further results are of interest, and are capable of straight- 
forward interpretation: 

(A) When (say) b becomes very small, there is no "reduction" in Qz, 
which is ~ 1 except in a small range of the strengths. In the limit of 6 = 
one has the simple polaron case, with no off-diagonal interaction. 
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Reduction Factors 
1 




Figure 1: Reduction factors vs. dominant coupling strength a. The diagonal 
Qz and the two ofF-diagonal reduction factors Qx and p are shown. Unlike 
the cubic-symmetry case, when q^ = Qz, in the dihedral case shown here the 
two Q'- factors are dissimilar, while p remains qualitatively similar to that in 
cubic symmetry. 
The ratio of coupling strengths - is maintained a constant, 0.8 



Reduction Factors 
1 




10 12 14 



Figure 2: Reduction factors vs. dominant coupling strength a. Same as 
figure 1, but with the ratio of coupling strengths - decreased to 0.667. The 
dissimilarity between Qz and Qx increases. 
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(B) When (instead of the cases shown in the figure) b > a, the roles of 
the two reduction factors Qx and qz are reversed. 

This can be justified as fohows: Interchanging a and b in the Hamiltonian 
of equation (4) has the effect of interchanging the operators ax and u^ . This 
interconverts, per definition in equation (23) and in equation (24) , Qx and 
Qz- Formally, the interchange can be performed by applying the unitary 
transformation matrix -7^{<^x + o'z) on the Hamiltonian. 

(C) The second relation in equation (27) , which is expected to hold 
for linear coupling to a single-mode coordinate in O symmetry, does not in 
general hold for neither q separately, but holds accurately for their mean. 
Thus we find under linear 6i , 62 coupling in dihedral symmetries that our 
data satisfy the new relations 

p = qx + qz-l (28) 

(D) Because the vibrational quanta are absent from the (proposed, 
quasi-classical) vibronic wave- functions, they do not play a role in the re- 
duction factors, as long as the coupling constants are defined in the non- 
dimensional form, in the way done here. 

4.2.1 Expressions for the reduction factors in dihedral symmetry 

The following expression relates the reduction factors to the H^-coefficients 

(29) 




defined and listed in [24] for point groups. (The definition of W allows 
various rearrangements of the symbols, not detailed now.) 
For electronic states belonging to a doublet E 

K{a) = (-l)WA(i?) Y.i-^PW {lll]<^.> (30) 



h 



b E E 



In the octahedral or dihedral groups ( — l)*-") is —1 for a = A2 and 1 oth- 
erwise. A(-E') = 2 is the dimension of the E representation. < x^ > is the 
weight of the nuclear component having the ^-representation in the vibronic 
state. The above formula is to be compared to the expression in Eq. (4.7.5) 
shown in [4], whose derivation is based on the Wigner-Eckart theorem, and 
to the corresponding expression for triplets in [31]. 
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Using equation (30) and Table D3.3 in [24] one arrives at the following 
expressions for the reduction factor in dihedral symmetry. 

K{Ai) = 1=<xa,> + <Xm> + <xI,> + < xl, > (31) 

KiA2) = p =< xm> + < Xm> - < xl,> - < xl, > (32) 

K{B2) = q, =< Xm>-<X%>-<xI,> + < xl, > (33) 

K{Bi) = q^ =< Xm>-<Xa2> + <xI,>-< xl, > (34) 

K{E) = (35) 

One notices immediately that K{B2) = qz i^ qx = K{Bi), as shown by the 
computed results in Figures 1 and 2. Similarly, upon adding up the first two 
lines and the last two lines separately and subtracting the sums from each 
other, one obtains 

p-iqx + qz-l)=^<X%> (36) 

This is similar to Eq 4.7.14 in [4] obtained in O symmetry. The right hand 
member is non-negative. However, for what are termed in [4] "ideal cases", 
the right hand side is zero and one recaptures equation (28) , as also found 
in our computation. Non-ideal cases are systems with coupling to more than 
one pair of modes [32], and others [33, 34]. 

4.3 Diagonalization within the vibronic doublet 

Two cases are of interest here. First, when the vibrational energies of the 
modes in equation (4) are finite (this excludes acoustic modes in a solid), 
and the external fields are weaker than the vibrational energies |il| << u;„ 
(all n). Then the admixture by the external fields of higher vibronic states 
can be neglected and one can work within the ground state vibronic doublet. 
This is carried out here. 

Secondly, when the external field components are periodic with a period 
-^ that is long, in the sense oftoj << |il|, so that during a period the system 
will stay in the lowest vibronic state (the adiabatic theorem) [17, 18]. Then, 
again, one can work solely within a ground vibronic state . This case will 
be treated in a future work. 

In terms of the reduction factors, the external field Hamiltonian in equa- 
tion (5) changes, as follows: 



n^ 


h^ 




-n-a- 


■^-^v 


• cJv, 


2 


2 





Hel-f = — ^ • I? -^ —^v ■ Ov = H^-f (37) 
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where 

Qv = iqx^x,P ^y, Qz^z) (38) 

and d\, are 2x2 Pauli matrices defined in the function space of the vibronic 
doublet. 

The constant interaction Hamiltonian H^^f has to be diagonahzed be- 
tween degenerate eigenstates of H^ in equation (7) . The diagonahzation 
sphts the states by an amount of 

2A = 2n\ny\ = 2hRy (39) 

and selects the following two linear combinations of the 9 and e states: 

|/ > = -sm — e 2 1^} + cos — e 2 |e} (40) 

\u> = cos—e 2 1^} + sm— e* 2 |e} (41) 

where we have defined two angles involving the three vibronic reduction 
factors and the three components of the field, through the expressions 

«. = „„,„„feM±(>Liwai ,42) 

and 

(/)v = arctan — -^ (43) 

Qx ^^x 

(The reason for the notation is that {R^,6^,(pv) make up the spherical co- 
ordinate representation of the "reduced" field vector Q^ in equation (38) 
.) |/ > and \u > are respectively the lower and upper split states of the 
doublet, but for $7^ < one has to take the branch between ^ and tt in the 
inverse tangent Oy. 

4.3.1 A counterintuitive effect of the off-diagonal coupling on the 
energy splitting 

For any given strength of the diagonal coupling, as the off-diagonal coupling 
increases, Qz decreases. This is seen by comparing in figure 3 the three 
curves depicting qz, (computed with a > 6, so that the diagonal coupling is 
dominant) in which the curves decrease as the off-diagonal coupling strength 
increases. This behavior is unexpected for the following reason: 
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q. 




Figure 3: Dependence of the diagonal reduction factor q^ on the off-diagonal 
coupling strength h. Reading the curves (here plotted against the dimension- 
less diagonal coupling strength a) from above to below, these were plotted 
for 6(< o) = 0.2,0.5,0.8 

When Q.X = ^z = 0, the splitting of the vibronic doublet components 
is, as we have just seen in equation (39) , 2|g2il2|. Here Vtz is a constant, 
while the multiplier qz is a function of the coupling strength parameters. 
However, the often quoted phenomenon known as the "repulsion of neigh- 
boring energy levels by interaction between them" would seem to require 
that qz should grow as the off-diagonal coupling increases, so as to make the 
splitting wider. This does not happen, and the reason is that the off-diagonal 
coupling hnq-n is not just a constant (a "magnitude"), but has a dynamical 
character. (We have purposely chosen for h low values, to show that the 
result shown here is not a high-order effect in h. We also recall some related 
discussion in the literature that lay the claim that the "repulsion between 
levels" need to happen only when the set of states is complete.) 

5 Application to Some Low Dimensional Systems 

5.1 Introductory remarks 

A number of organic hydrocarbon molecules or radical ions provide instances 
for a doublet interacting with two independent modes. Investigations of their 
stereochemical properties (that is, their possible configurations) in either 
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the ground or excited states have been lucidly summarized by Bersuker [35] , 
giving also extensive references. 

When the molecules belong to a dihedral symmetry group, they can be 
subject to linear couplings by two types of non-totally symmetric nuclear 
displacement modes. When these are of unequal strength, the adiabatic po- 
tential surface will have two minima along the dominant coordinate and two 
saddle points along the non-dominant one. [Couplings quadratic or of higher 
order in the mode coordinates (not considered in this paper and also not of 
paramount importance in many hydrocarbons) can turn the saddle points 
into minima [30].] Further distortions from these simple configurations are 
also possible [35, 36]. Typical energy differences between alternative con- 
figurations are of the order of 2.5 kilo-kelvins, which are also the values in 
these systems for the Jahn- Teller (or stabilization) energies that are relevant 
to the considerations in this paper. 

Transitions between minima take place typically across the saddle point 
(a "transition state"); this is discussed for the cyclobutane radical cation 
IC4H-+) in [37]. (Figure 4.) 

The activation energies for the transitions are estimated in the same 
source as perhaps 1 kilo-kelvin, possibly dropping below one half kilo-kelvin. 
With dynamic processes included, as in the present paper, the "transition" 
is of course an ingredient of the ground state wave function, (also termed 
"tunnelling") and does not represent a real process. Exceptions are when 
the molecule is embedded in a matrix with lower symmetry than the nominal 
one. Then the molecule may be forced into one of the minima, from which 
it can make a real transition to another minima. A further instance of real 
transition is a thermally activated one. The two states obtained in the last 
section, [equation (40) , equation (41) ] can be the basis for calculating 
transition rates but, typically, in a thermally activated process one includes 
states higher than the ground doublet and their inclusion is outside our 
concern here. 

Molecular data (like potential energy surfaces) have been calculated for 
the cyclobutane radical cation in [37], for the cyclobutadiene {C4H4) radical 
cation in [36] , and for radical cations of several cycloalkanes in [38] . We cite 
[35] for other systems. Calculated results tend to be very sensitive to the 
level of the computational effort (and this sensitivity also includes the order 
of relative stability among different configurations). In particular, single 
determinant wave-functions appear to be unreliable. 

Electron spin resonance (ESR) experiments have the capability of throw- 
ing light on dynamic effects, which we have studied here. Experimental de- 
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D2d-D2 








b2 

-4- 











Figure 4: Distortions and states of the cyclobutane radical cation (CiH'g^). 
The upper part shows two possible distorted configurations of the frame of 
the four carbon atoms. Experiments favor the B configuration {C2v)- The 
puckered form of this configuration(two carbon atoms above the paper and 
two below) is not shown; neither are the positions of the eight protons (also 
outside the plane of the paper). The lower part shows the splitting of the 
energy levels of the highest lying partly occupied doubly degenerate e-state 
after distortion, as well as the shapes of the orbitals. Three spinning elec- 
trons are placed in the levels. The three-electron, product state is labelled in 
the text 9, while the state in which 62 is doubly occupied and 61 singly occu- 
pied is labelled e. Note that the electron densities on two-two carbon atoms 
differ in B, but are the same in A. (With permission, after [39]. Copyright 
(1983) American Chemical Society.) 
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termination of the vibronic reduction factors (qxiPiQz) from the observed 
5f-factors in the spectra are unfortunately unhkely, due to covalency effects 
and the smah spin-orbit couphng in these cychc compounds ( this is unhke 
the transition ion compounds [12, 28]). However, the behavior of the proton 
hyperfine hues in hydrocarbons can give a clue to dynamic processes and, 
especially, to their coalescence and narrowing as the temperature is raised 
[38, 39]. 

5.2 ESR in the cyclobutane radical cation 

In particular, we wish to apply the present theory to the hyperfine lines of 
C^H'^ in a solid matrix observed and discussed by the authors of [39]. Their 
measurement of the electron spin resonance absorption at and above 77^ K 
gave fairly equidistant nine-line spectra, which are indicative of electron- 
nuclear interaction with eight equivalent protons. On the other hand, low 
temperature observations at A.2^K showed hyperfine lines with different 
separations. The observed separations are consistent with non-equivalent 
coupling to the four pairs of protons. This non-equivalence can be under- 
stood if one has a puckered molecule in which the four carbon atoms have 
undergone further distortion from D2d to C2v This entails a 62 distortion- 
mode leading to a rhombus-like structure, as illustrated in figure 4. In our 
terminology, for the coupling coefficients this implies that \a\ > \b\. 

Part of the non-equivalence among the proton pairs may well disappear 
at higher temperature, either by the flattening of the puckered molecule 
or by thermal averaging between alternative non-planar forms (assuming a 
softening of the barrier in the matrix at the higher temperatures). However, 
the remaining distinctness between two pairs of carbons (located respectively 
on the blunt and the sharp angle apexes of the rhombus) cannot be fitted to 
the high-temperature isotropic spectra. It has been suggested in [39], that 
the dynamic Jahn- Teller effect is capable of removing this non-equivalence. 

Since the quasi-classical theory enables one to quantitatively evaluate 
the extent of "democratization" of the carbon atoms, we shall now provide 
some numerical estimates for this. 

5.2.1 Electron densities 

The effects of Jahn- Teller distortions and of the motions between them on 
the hyperfine structure have been clearly formulated in [40, 41] and amply 
reviewed in [42], so that we do not have to repeat the theory. In essence. 
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what one sees is the proton nuclear spin interacting with the electronic spin 
density on the carbons. For either of the carbon pairs, Ci and C3 on the 
sharp angled apexes and C2 and C4 on the blunt angled apexes, the spin 
density will differ in the Qe and Ce electronic states in equation (1) . The 
density-difference of the two states will express itself in the spectrum in an 
opposite manner. 

On the other hand, the nuclear function cofactor of C,e and Ce in the 
vibronic state will regulate the relative weights of these electronic states in 
the vibronic state. We are now looking for the average weights of these 
electronic states. They are clearly given by the expectation values of the 
projection operators Pq and P^ in the vibronic state. These electronic pro- 
jection operators can be written as 

P,/, = i(/±a,) (44) 

Suppose now that the host matrix of the cation radical stabilizes the Ce 
electronic state in preference to the C,^ state. This stabilization can come 
about by having the following external field parameters : O^, = 0^ = and 

JI2 = -|$7°| < 0, leading to 

H,i^f = -^\QP\a, (45) 

(Cf. equation (5) .) In the absence of any coupling to the nuclear motion 
the pure electronic C,^ state is stabilized. This state leads clearly to non- 
equivalent couplings, whereas equivalence is regained only for states in which 
Ce and Ce appear with equal weight. 

Let us now turn on the coupling to the nuclear coordinates. By the 
results in section 4.2, the lower energy solution is equation (40) . Since, 
with our choice of the external fields, the angle 6^ = tt, this solution is 
simply the vibronic state |^} shown in equation (15) . We can now evaluate 
the expectation values of the electron projector operators in this state by 
using the definitions of the reduction factor qz in equation (24) . We get 

{9\Pe\e} = 1(1 + qz) and {9\P,\0} = ^(1 - g,) (46) 

What happens in the dynamic case? Remarkably, even here one does 
not achieve the extent of uniformity in the spin density, which is achieved in 
the high symmetry case (e.g., O) for which q^ = 0.5, since, as seen in figures 
1 and 2, q^ decreases only slightly below 1, so that the difference between 
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the two projectors will still remain. (Of course, this comes about because of 
the dominance of the (rhombic) 62 distortion over the vying 61 distortional 
mode.) The end-result is that under realistic conditions that the (rhombus- 
inducing) 62 coupling is significantly stronger than the (rectangle-making) 
61 mode, the DJTE cannot make the two pairs of carbons equivalent, or the 
spectrum equidistant. For this to happen (at some elevated temperature), 
either the external splitting field Q^ must be lowered, or fast jumps between 
the two vibronic eigenstates |^} and |e} or to higher lying vibronic states 
will have to occur. ("Fast" means shorter than 10~^ seconds, the hyperfine 
coupling time scale.) 

5.2.2 Dynamic effects on hyperfine lines 

A numerical estimate for the c/z-factor confirms this conclusion. We have 
first estimated the coupling strengths for two coupling modes in the cyclobu- 
tane radical cation from computed stabilization energies in the rectangular 
and rhombic configurations, shown in Figure 7 of [37]. The values obtained 
from the UMP2/6-31G*//UHF/6-31G* variational method were adopted, 
since this places the rhombic configuration below the rectangular one, as is 
observed. The computed stabilization energies relative to the square con- 
figurations are 6600 cm~^ (for the rectangular shape) and 7250 cm~^ (for 
the rhombic form). Taking the experimental wavenumbers in cyclobutane 
observed by [43] for the modes: namely, 926 cm"^ in the C — C stretching 
(61) and 1001 cm~^ in the CCC angle bending (62) modes, we obtain for 
the dimensionless coupling strengths: 

Ob2 = 7.611, h^ = 7.551 (47) 

We note that these values are close to each other. 

From our expression for the vibronic wave function and a simple quadra- 
ture, we obtain for these values a reduction factor 

q^ = 0.999998 (48) 

that is, near enough unity, thus excluding the attainment of equidistant 
spectra by a pure dynamic mechanism alone. Phrased alternatively, the 
relevant matrix elements of the relaxation matrix are proportional to the 
overlaps between Gaussian wave functions localized around different minima. 
These are greatly reduced by the strong vibronic coupling [44]. 

We recall that in O symmetry, when the two coupling strengths \a\ and 
\h\ are precisely identical, so that Qz = q-x = Q, for strong linear coupling these 
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reduction factors approach 2, which indeed halves the difference between 9 
and e occupancies. In D2d symmetry, when the couphng strengths a and 
b are different and large >> 1, the tunnelling between different stabilized 
wells is negligible. 

6 Suminary and Outlook 

We have tackled (fully, though not exactly) the time-independent quan- 
tum mechanics of a pair of isolated electronic states, and this under rather 
general conditions: namely, the states are subject to interaction both with 
static external fields and with a dynamic surrounding, the latter in the lin- 
ear approximation. Our non-perturbational treatment was made possible 
(a) by having a closed solution (the quasi-classical vibronic wave-function) 
for the part expressing the coupling between the electron and its dynam- 
ical surrounding, and (b) by inverting the usual order of solution through 
taking step (a) first and including the external field later. Thereby, the 
dynamically coupled states maintain convenient symmetry-group properties 
in the (Hilbert) function space; the external forces are subsequently treated 
within this framework. Their strength is, however, modified ("renormal- 
ized") by "vibronic" reduction factors. The use of these factors in a non- 
perturbational way is yet another new feature of this approach. 

Branching out from the present treatment carried through for two states, 
one can similarly tackle problems with interactions affecting three arbitrary 
electronic states, as well as any chosen number of states. By extension of 
the approach worked out in this work and shown in equation (11) , one 
finds in these cases also that the wave-function generator, the generalization 
of ^, is of a closed form. This consists of a finite number of terms, each 
with a given symmetry. Precisely, the generators (matrices) correspond to 
all representations of the symmetric product of the electronic multiplet in 
its reference group. Thus, recalling the results in this paper, for a doublet 
we have three matrices (those appearing in equation (13) ). (These three 
matrices are recognized as representatives in O of the symmetric product 
of the E-representation, namely Ai and E.) In a similar manner, for an 
electronic triplet one has six matrices: namely, the identity matrix, three 
angular momentum matrices and two more trace- less matrices. All these 
make up the generator wave-function. Likewise, for any larger number of 
states. Algebraic relations connect the functions belonging to each term, 
which relations have to be solved simultaneously. (The point of these re- 
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marks is to assert that the doublet is not a fluke-case, but rather a special, 
though by far the simplest, case of multiple electronic states in interaction 
with their surroundings.) 

In conclusion, we restate that the major restrictions on the applicability 
of this work and of its possible extensions are the validity of regarding a 
finite number of states in isolation and the linear approximation. 
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